Spin relaxation in isotopically purified silicon quantum dots 



M. Prada 1 ' 2 , R. H. Blick 2 , and R. Joynt 1 

1 Department of Physics, 1150 University Ave., University of Wisconsin-Madison, Wisconsin 53706 (USA) and 
^Electrical and Computer Engineering, University of Wisconsin-Madison, Wisconsin 53706 (USA) 

We investigate spin-flip processes of Si quantum dots due to spin-orbit coupling. We utilize the spin-orbit 
coupling constants related to bulk and structure inversion asymmetry obtained numerically for two dimensional 
heterostructures. We find that the spin-flip rate is very sensitive to these coupling constants. We investigate 
the nuclei-mediated spin-flip process and find the level of 29 Si concentration for which this mechanism become 
dominant. 



I. INTRODUCTION 

The control of single-electron spin in nanoscale devices is 
a key element in the field of spintronics, where the spin de- 
gree of freedom is used for information transfer and process- 
ing. The vast majority ofproposed spin devices are based 
on semiconductors |Uj,|2l|3t|4j], where single-electron spin has 
been probed to be controllable by means of voltages applied to 
electrostatic gates [5]. A promising technology for the imple- 
mentation of quantum computation (QC) involves the storage 
of quantum information in the spin of electrons in quantum 
dots (QDs). The key requirement is that the lifetime of the 
spins is long compared with the time required for the oper- 
ation of logic gates. This has motivated the development 
of QDs in Si-based materials [6], where spin-orbit coupling 
(SOC) is weak and isotopic enrichment can eliminate hyper- 
fine coupling (HC). In such scenario, processes limiting QC 
are dominated by SOC resulting from spatial inversion asym- 
metry. In typical solid state systems, macroscopic electric 
fields cause structure inversion asymmetry (SIA), giving rise 
to Rashba-type terms [8] of the form: 

H R = a(6 x k y -a y k x ), 

whereas fields resulting from the lack of an inversion center 
lead to bulk inversion asymmetry (BIA) and to the Dressel- 
haus term, H : 

Here, o, and ki denote spin and momentum operators, respec- 
tivelly, and a, p are the Rashba and Dresselhaus coupling con- 
stants, respectivelly. Although both SOC contributions have 
been noted for decades, their absolute value have been mea- 
sured simultaneously only recently [11]. The different the- 
oretical models show controversy [12], with debated estima- 
tions for Rashba or Dresselhaus related parameters. 

In this paper we calculate spin-flip rates for typical Si QD 
via higher-energy virtual state, involving also phonon emis- 
sion. We define the regime for which SOC is the dominant 
source of relaxation as a function of 29 Si isotope concentra- 
tion. We utilize in our calculations the SOC parameters ex- 
tracted numerically using for first time a model that includes 
atomic crystal symmetry, spin and interfaces built into the ba- 
sis representation 01311 . that aims to clarify the existing con- 
troversy. 



This paper is organized as follows: Sec. ITTI contains a de- 
scription of the model: we describe the electron-phonon in- 
teraction, and the spin-flip mechanism, namely SOC and hy- 
perfine coupling (HC). In Sec. |III]we present our results, and 
finally, Sec. [IV]is devoted to conclusions. 



II. METHODS 
A. Model 

We consider a QD formed in a two dimensional (2D) het- 
erostructure with parabolic confinement. Spin-flip is consid- 
ered via an orbital (or valley) state with energy fixoo, as de- 
picted in Fig. Q] 
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FIG. 1 : Schematic representation of a spin-flip process for a Si QD with one 
electron. 



In order to calculate the relevant transition rates for spin 
relaxation, we consider a perturbation given in general by 
&H = H s f + Hph, 7/ph corresponding to the electron-phonon 
coupling, and H s { to the mechanism causing a spin flip in the 
dot. The amplitude of a spin-flip event in the ground state 
then is given in second order perturbation theory, 
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where A is the energy exchanged with the bath. Here, 1, 1 de- 
note the spin state and n is the orbital degree of freedom. The 
orbitals are described in terms of Fock-Darwin states. The 
ground-state spin-flip rate to |1 1) is then given by Fermi's 
golden rule: 
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In this notation, |1 J.) denotes the initial state of electron, nu- 
clei and phonons, |1 I) — » \i) = |1 J.) ® \i n ) ® |/ pn ), likewise 
|1 T) -> I/) = |1 T) ® |/») ® l/ph). Inserting Eq. © into © 
we obtain (setting A — >0): 

?u = xl<^I^Ml 2 «(^-g/) K1T ff i>|2 - (3) 

We note that the phonon and the spin-flip parts can be factor- 
ized, so we define Ft i = r pn x i^ 2 , with 
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Hence, r pn and £ can be treated separately. r p h describes 
the phonon-emission process, coupling two orbitals via 
electron-phonon deformation potential, and £ includes the 
spin-flip process. We give next a detailed description of both 
terms. 



B. Phonon-electron coupling 

For Si under compressive stress along [001], the electron 
interacts with a phonon of momentum q via deformation po- 
tentials HEUE]- The Hamiltonian reads: 
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[a*e iqr + a qs e iqr ]q{Z d ^ x q x + Z d e%, - 
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(Z d + Z u )e s z q z ) 



(5) 



where (n q — \ \a q \n q ) = (hn q /2M c (O q ), M c is the mass of 
the unit cell and h& q the phonon energy. Here, s denotes the 
polarization of the phonon (two transverse and one longitudi- 
nal), q is the wavevector, Z u and Z d are the electron-phonon 
coupling parameters in Si (see table Q). This is slightly sim- 
pler than the corresponding Hamiltonian in GaAs because of 
the absence of the piezoelectric coupling in Si. In order to 
calculate r p h, we use the electric dipole (ED) approximation, 
e iqr fa 1 + iqr (valid in the range of 0.1-0.2 meV, which is 
roughly the phonon energy involved in this relaxation pro- 
cess). 

For the process depicted in Fig. Q] mixing of s-orbitals occur 
at 0-th order in ED (to lowest order). We summarize in ta- 
ble|T]the parameters relevant for the computation of TV i used 
throughout the text. Using the 0-th order ED approximation, 



S„[eV] 


E d [eV] 


p si [kg/m 3 ] 


v,[m/s] 


v,[m/s] 


A [/jeV-nm 3 ] 


<p 2 )[nm 2 ] 


9.29 


-10.7 


2330 


9330 


5420 
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TABLE I: Physical constants and material parameters for Si and the QD as 
in 01. 

e iqr « 1, we obtain for the electron-phonon coupling: 
#ph ~ E( fl ? +a\) W\ [Zde x q x + Z d e y q y + (Z d + Z u )e z q z ] , 



which, using ©, gives the rate: 
K + l) , 
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Here, n q +\ « 1 in the range of low temperatures considered, 
and psi is the density of Si (see table e, is the 2-th com- 
ponent of the unitary polarization vectors and s denotes the 
polarization, s=\, ti, t2'. 
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and q = (sin8cos(p,sin8sin(p,cos8). To integrate Eq. ©, 
we assume an isotropic phonon spectrum, = HG> qs with a 
dispersion relation <x> qs = v s q, v s being the sound velocity of 
the mode s. This gives: 
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Using the above result, the polarization vectors of (0), rear- 
ranging (O, and performing the angular integral, we obtain a 
compact expression for the rate: 
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C. Spin orbit coupling 

As pointed out before, spin-flips can be provided by two 
different mechanisms, the HC with the 29 Si nuclei and the SO 
coupling. For an isotopically-purified sample, we expect SO 
to be the dominant spin-flip mechanism. Hence, we first eval- 
uate SO mediated spin-flip and then consider the 29 Si concen- 
tration at which the HC mediated spin-flip rate becomes com- 
parable to the SO mediated one. 

We express the Hamiltonian in the convenient phase coor- 
dinates (,quqi,puP2) B: 
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with C0i,2 = (cOq + CO 2 ) x l 2 ± C0 C and co t = eB/m. 

We add as a perturbation the Rashba and the Dresselhaus 
SOC terms in a convenient basis rotated by 45° (see Fig. |2j: 

#so = ~ [(P + a)d+/5_ + (p - a)o_/3+] . (10) 

We note that the SOC term does not commute with Hq. This 
implies that the spin does not remain fixed during the motion 
of the electron, apart for an apparent homogeneous shift of the 
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FIG. 2: (a) Schematic representation of a spin-flip process for a Si QD with 
one electron, (b) Spin direction for the eigenfunctions of Rashba and Dres- 
selhaus contributions: for p || [110], [110], the effective magnetic fields are 
parallel. 



momentum to first order, which would not change observables 
(other than just a shift in the energy levels). To get rid of 
such terms fixed by gauge invariance, it is convenient to apply 
the transformation of Aleiner et al, [ 17] U = exp{/yi^ + d_ + 
iy 2 q-6+}, (q± is the position operator along the directions 
e±), which allows to see how the levels split in the presence of 



Hso- Using the Zassenhaus formula, e A+B : 
we factorize U to second order in SOC: 
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Retaining terms to third order in SOC, the momentum trans- 
form as: 

U^P+U ~ P+ + hjid- - hjiy 2 q-O z - 2kfy 2 q-q+d+, 
and likewise: 

U^P-U ~ P + hy 2 a + + hy^q+Oz - ihy^q-q+o- . 

so that the total Hamiltonian, H = Hq + Hso, transforms as 
H' = WHU. Considering terms up to third-order in SOC: 
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where we have defined y\=m(a — P)/ft = l\ and y 2 = 
— m(a + P)/fi = —V}) with h. 2 characterizing the length 
scale associated with the strength of the SOC for electrons 
moving along the crystallographic directions. This choice al- 
lows us to cancel out the linear terms in momentum. The 
second term of Eq. (fT~2T > gives an effective magnetic field that 
can be expressed as: 

H\ = ih\\a\al — a y d\\o z , (13) 

where we have defined: 

V V <*>! V 6)2/ 

We note that Hi does not break the Kramers degeneracy, and 
hence, does not contribute to spin-flip. However, the last term 



of Eq. ( [T2l i does break Kramers degeneracy, allowing spin-flip 
terms. In second quantization, we have: 
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We then conclude that, for SOC, i^so ~ if / ^CO, from which we 
can proceed to calculate T| |, using < fT~4T > and ((HJ in (0J. We 
emphasize that this term is third order in SOC, and applying 
Fermi Golden's rule, the spin-flip rate will appear in sixth or- 
der. Hence, T-t 1 is very sensitive to the SOC parameters, so 
it is critical to determine them to a high degree of precision. 



D. Hyperflne Coupling. 

For a electronic spin in a QD and in the presence of nuclear 
spins, a contact interaction Hamiltonian can be assumed as a 
perturbation, 



H^^^MBMmSiljSirt-Rj) = A^5,/ ; -5(r, -Rj) 

i.j U 

(15) 

where 5, {If) and r, (Rj) denote the spin and position of the / th 
electron (/ h nuclei), and A is the hyperfine coupling constant. 
We now use ( fT3] > as the spin-flip mechanism to calculate ^ of 
Eq. ©: 

Kltltfhcl^l^A^KT \S + IT8(r-Rj)\l)\ 2 . 



Next we substitute the sum over j by an integral over the nu- 
clei with a density C„: 

KiTl^e^oi^A^ld^-io,^)! 2 ^^)! 2 

which can be easily performed assuming Fock-Darwin states 
for <!>;, resulting: 
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Using ( fTSI l together with Eq. dH) in (0), we can thus obtain the 
spin-flip rate due to HC. Typical values for r| 5 ' of 186 have 
been reported |H d yielding A « 2 x l(T 7 eV-nm 3 . 
Only about 4% of the nuclei have spin, so C n « 0.04 x 8 /vo « 
2nm~ 3 (vo » 0.17nm 3 , is the unit cell volume for Si). It is im- 
portant to note that is proportional to C„, i.e., to the total 
number of nuclei N„ with which the electrons interact. This 
is consistent with the simple picture that the relaxation rate 
is proportional to the mean-square fluctuations in the random 
hyperfine field. Formulas for spin relaxation rates due to HC 

-1/2 

that give an apparent proportionality to N n are common in 
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the literature, and have given rise to the incorrect notion that 
some sort of motional narrowing is at work. This is not possi- 
ble, since the fluctuations in the nuclear spin system are slow. 
In any case the rate must vanish as N n — > 0. These formulas 
are correct, but they generally involve other parameters that 
actually vary with N n . 



responds to the horizontal solid line of Fig. [3] We can see that 
a purification of about 80% of the 29 Si isotope will cause SOC 
to become dominant. In contrast, for a similar QD formed in 
a Si membrane, with purification of ~ 99% we obtain SOC 
as the dominant mechanism, increasing the relaxation time by 
almost two orders of magnitude. 



HI. RESULTS 

We considered a typical QD formed as in |01 with a level 
separation of h&o = 0.2meV. The a and p parameters where 
extracted numerically using NEMO-3D B26I1 on nanoHUB.org 
computational resources B27I1 . In NEMO-3D, atoms are repre- 
sented explicitly in the sp 3 d 5 s* tight-binding model, and the 
valence force field (VFF) method is employed to minimize 
strain [28]. NEMO-3D enables the calculation of localized 
states on a QW and their in-plane dispersion relation with a 
very high degree of precision, allowing to extract the split- 
tings along the in-plane directions in momentum space 11311 . 
For now, we note that the a value is well defined and depends 
on external electric fields applied to the sample, whereas the 
Dresselhaus, (3, depends strongly on the atomistic details of 
the interface. 
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FIG. 3: Tji as a function of 29 Si isotope purification degree. The hori- 
zontal lines correspond to T-t i due to SOC, for various values of Dresselhaus 
coupling parameter, p. 

Fig. |3]shows the calculated values for Tj j via both mecha- 
nism, HC and SOC. The parallel lines correspond to r^j 50 for 
different possible values of p, whereas the thick solid line rep- 
resents the spin-flip due to HC mechanism, r^ c , as a function 
of 29 Si isotope purification degree. Using our numerical value 
for a and a few different values for p, we find that for natural 
Si, the HC dominates. We note that these depend strongly on 
the sample, so a Si sample with dominating SOC is possible. 
In particular, for a Si/SiGe heterostructure, the r^ c value cor- 



FIG. 4: Regions where SOC (shadded) or HC (light) are the dominant mech- 
anism of relaxation as a function of Dresselhaus parameter, (5, and degree of 
29 Si isotope purification: for small p, a larger isotopic purification degree is 
needed, and hence, larger relaxation times can be achieved. 

Fig. |4]represents the regions for which SOC or HC are the 
dominant spin-flip mechanism. It also summarizes the most 
relevant P values found numerically: a typical QD defined in a 
Si/SiGe heterostructure has p ~75m/s, for which HC is dom- 
inant up to a purification level of 81.5%. A similar structure 
obtained on a Si membrane will show a larger Dresselhaus, 
and hence, larger purification will be needed to have SOC as 
dominating mechanism. 



IV. CONCLUSIONS 

In summary, we have calculated the dominant rates for 
phonon-assisted spin-flip relaxation for the ground state of a 
single-electron Si QD. Tf j is found to be of the order of tens 
of seconds, very sensitive to SOC. We observe that choosing 
a sample with small SOC, as pure Si membranes, can lead to 
spin lifetimes of up to a few hours. QDs fabricated in Si/SiGe 
heterostructures are more limited by SOC, and purifications 
of only about 80% would give spin lifetimes of a few minutes. 
Due to weak spin-orbit and HC, Si offers very long coherence 
times, which are required for solid state qubits. 
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